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Abstract 

A simple mathematical procedure is introduced which allows redefining in an 
exact way divergent integrals and limits that appear in the basic equations of 
classical electrodynamics with point charges. In this way all divergences are at 
once removed without affecting the locality and the relativistic covariance of 
the theory, and with no need for mass renormalization. The procedure is first 
used to obtain a finite expression for the electromagnetic energy-momentum 
of the system. We show that the relativistic Lorentz-Dirac equation can be 
deduced from the conservation of this electromagnetic energy-momentum plus 
the usual mechanical term. Then we derive a finite lagrangian, which depends 
on the particle variables and on the actual electromagnetic potentials at a 
given time. From this lagrangian the equations of motion of both particles 
and fields can be derived via Hamilton's variational principle. The hamil- 
tonian formulation of the theory can be obtained in a straightforward way. 
This leads to an interesting comparison between the resulting divergence-free 
expression of the hamiltonian functional and the standard renormalization 
rules for perturbative quantum electrodynamics. 

03.50.De, ll.lO.Ef 



Typeset using REVTeX 



*Email address: massimo.marinoOmi. infm.it 



1 



I. INTRODUCTION 



The relativistic Lorentz-Dirac equation was derived in 1938 [|I| to describe the motion 
of a classical point charge in an arbitrarily assigned external field, taking into account its 
electromagnetic self-interaction and therefore the effect of radiation reaction (useful reviews 
on this subject can be found in Ref. 

Dirac believed that, once a finite classical theory for point-like charges had been ob- 
tained, then by applying appropriate rules of canonical quantization it should be possible to 
formulate the quantum theory in a way that avoided the occurrence of troublesome infinities. 
Already in the final part of his original paper [|I|, with the explicit purpose of suggesting 
the path towards a possible quantization, Dirac showed that his theory could be deduced 
from a variational principle, by making use of a non-local action which had been previously 
proposed by Fokker [|] to describe a system of charged particles interacting at a distance 
through the average of the retarded and advanced fields. However Dirac soon realized that 
such an action could not allow to reexpress the djTiamics in hamiltonian form. For this 
reason, soon thereafter he produced a second paper ^ specifically devoted to the problem 
of achieving a canonical quantization of his theory. In this work he abandoned the action 
principle, and attempted to write down directly a four-dimensional hamiltonian in terms of 
an unphysical field dependent on the proper time of each particle constituting the system 
under consideration. For this field Dirac skilfully postulated particular ad hoc Poisson brack- 
ets, involving an arbitrary time-like constant four- vector A, in such a way that his covariant 
equation of motion for the particles was recovered in the limit of vanishing A. It is not 
surprising that such an involute approach could hardly lend itself to further generalizations, 
and Dirac was unable to use it as the starting point for a consistent formulation of quantum 
electrodynamics. It is also interesting to observe that, although Dirac had derived in his 
equation of motion using essentially the principle of local energy-momentum conservation, 
he did not provide a general expression for the global conserved quantities. Besides, the 
conservation laws no longer emerged in an obvious way from the canonical formalism that 
he proposed. 

With the advent of renormalization theory, perturbative quantum electrodynamics devel- 
oped in an impressive way without relying on any input from classical theory for the formal 
treatment of divergences. However, research on the classical theory was never abandoned. 
In particular, a new finite classical action describing the behavior of both particles and fields 
was obtained by Rohrlich by adding an appropriate term to the original expression of 
Fokker and Dirac. However, this result could not help in providing a canonical formulation 
of the theory. Furthermore, in the action the electromagnetic field was split as the sum 
of a source-free term plus the contribution generated by the particles through the retarded 
and advanced potentials. This contrasted with the lagrangian formulation of QED which, 
despite the subtleties related to the problem of providing an adequate renormalization, is 
always expressed only in terms of the full interacting electromagnetic field. Similar con- 
siderations apply also for other studies subsequently appeared, in which the expression of 
the electromagnetic energy-momentum in terms of the retarded fields was subjected to an 
accurate analysis 0, and a different kind of lagrangian, based on the distribution theory, 
was derived [Q. 

In the present paper we shall propose a simple and general method to evaluate the inte- 
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grals of the electromagnetic energy-momentum, lagrangian and hamiltonian densities over 
spatial regions containing the point particles. In a similar way, we shall assign a finite value to 
the fields and potentials on the particles' world-lines, by generalizing a procedure previously 
proposed by Teitelboim and Lopez In this way practically all the familiar formu- 

lae of electrodynamics, which are commonly applied when the particles' self-interaction is 
neglected, can be immediately generalized to the self-interacting case by appropriately rein- 
terpreting each individual term according to the given prescriptions. The resulting equations 
of motion will be the relativistic Lorentz-Dirac equation for the particles and the Maxwell 
equations with point-like sources for the fields. In particular we shall obtain: 

1. A finite expression for the total energy- momentum of the system of fields and particles, 
which is Lorentz-covariant and conserved on the solutions of the equations of motion; 

2. A finite lagrangian which determines the dynamics of both fields and particles accord- 
ing to Hamilton's variational principle; 

3. A finite hamiltonian and a system of Poisson's brackets which express the dynamics 
in canonical form. 

In the treatment here presented, the infinities will be consistently removed without alter- 
ing any of the fundamental formal properties of the theory, including locality and relativistic 
covariance. A remarkable property of this procedure is that it allows expressing the four- 
momentum, the lagrangian and the hamiltonian of the system only in terms of the particle 
variables and of the total instantaneous electromagnetic field, i.e. without any unphysical 
splitting into retarded, advanced, incoming or outcoming parts. This will provide a formula- 
tion of the theory which is in complete agreement with the requirements of local relativistic 
field-theory. Furthermore, no diverging parameter will have to be introduced, since only 
the physical mass of the particles will be used. Due to the close correspondence with the 
quantum theory, we shall show that these results may also suggest a new way of formulating 
the renormalization procedure of quantum electrodynamics. 

The paper is organized as follows. In Sec. II we define the renormalized energy- 
momentum associated with the electromagnetic field in the presence of point charges, and 
we discuss how the conservation of the sum of this energy-momentum and the familiar me- 
chanical term associated with the particles (with mass parameter equal to the observable 
mass m) determines the equation of motion of the particles, provided that the field satisfies 
the usual Maxwell equations. Actually, this provides a new rigorous and straightforward 
derivation of the relativistic Lorentz-Dirac equation. We shall comment on the mathemat- 
ical structure of the system of coupled equations of motion for particles and fields, and on 
the non-trivial constraints that are to be imposed on the initial data of the corresponding 
Cauchy problem. In Sec. Ill we define the renormalized lagrangian of the system, and we 
show how the equations of motion of both particles and fields can be derived from an action 
principle. We shall see that on the solutions of the field equations the value assumed by 
our lagrangian can be approximated by an expression corresponding to an electron of finite 
radius r and bare mass rriQ = m — q'^/Snr, although the canonical forms of the two theories 
remain essentially inequivalent even in the limit r — > 0. In Sec. IV we give a precise math- 
ematical definition for the partial derivatives of the lagrangian with respect to the particle 



3 



variables, and we make use of this definition to put the equations of motion in the stan- 
dard form of Euler-Lagrange equations. The reformulation of the theory in the hamiltonian 
formalism will then be obtained in a straightforward way via a Legendre transformation. 
In Sec. V we obtain an equivalent expression for the renormalized hamiltonian in terms of 
the Fourier transform of the electromagnetic potential, and in Sec. VI we make use of this 
expression to illustrate the close relationship with the renormalization rules for the electron 
self-energy in perturbative quantum electrodynamics. 

Most of the mathematical proofs and technical developments are worked out in the 
Appendixes. In particular, in App. A we consider the retarded and advanced fields near a 
point charge in arbitrary motion, and we show that they can be asymptotically expanded 
as the sum of a series of covariant homogeneous functions of the space-time variables. This 
essential result will be extensively applied throughout the paper. 



II. THE RENORMALIZED ENERGY-MOMENTUM 



A. The electromagnetic energy-momentum in the presence of point charges 

In the elementary case of electrostatics, the energy of the electric field associated with a 
non-singular spatially confined charge distribution p(x) is given by 

U = y£^E'i^), (2.1) 

where 

E(x) = / dVp(x')^^^. (2.2) 
An J |x — x'l 

When the charge distribution is constituted by a set of point charges qi, q2,- ■ ., Qn located 
at the points Zi, Z2,. . ., z^, Eq. ( |2.2| ) takes the form 

1 ^ x- z 

^" n=l |X Z„| 

It is well known that in such a case the integral on the right-hand side of Eq. ( |2.1| ) becomes 
divergent. However, introducing a cut-off at a distance r from the point charges, one can 
define 



— 1 r ^ 



N 

" ^'x - z„| - r) 

n=l 



where 6 is the usual step function. Then an elementary integration provides 

U^r) = ^^^:1l + ^T.T^^^+0{r) forr^O. (2.3) 

It can be seen therefore that the divergent part is only represented by an irrelevant constant 
term which can be safely dropped, leaving us with the familiar Coulomb potential energy: 
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t/(zi,...,z^) = i-E^^^- (2.4) 



n<m I'^n 



The above considerations illustrate the obvious fact that the singularities associated with 
point charges are completely harmless in the case of electrostatics, since the physics can be 
completely described by finite quantities such as (|2.4| ). In the present paper we shall show the 
remarkable fact that an analogous result actually holds also for relativistic electrodynamics. 
To this end, it is useful to introduce a general limiting procedure for functions near a first- 
order singular point. 

Definition 1 Given a function f{r) of a variable r, having in a neighborhood of r = a 
behavior of the form 

fir) = ^ + 6 + 0(r), 
r 

we define its "TZ-limit" (henceforth indicated by the symbol T^limj for r — > as 

TZ lim /(r) = lim \rf{r)] = b. 

r^O r— >0 dr 

Correspondingly, given a function G'(x) of the spatial coordinates x, having in a set V the 
singular points zi, Z2,. . ., zn, and a behavior such that 

r ^ a' 

/ d'^G{^) n ^ (|x - z„| - r) = - + 6' + 0(r), (2.5) 
we define its "TZ-integral" over V as 

N 

n d^xG(x) =niim d^xG(x)Y\ 0(\x-zJ-r) = b'. 
Jv '"^0 Jv 



It is clear that the 7^-integral defined above differs from the ordinary integral only in a 
neighborhood of the singular points. More precisely, if ^ = U V2, where Vi fl V2 = and 
V2 contains no singularities, then 

7^ / c/=^xG(x) = n f rf3xG(x)+ / d3xG(x). 
Jv JVi JV2 

Applying the Definition 1 to the electrostatic problem considered at the beginning, we 
can write 

U{z„...,z^) = ^Y. = nfd'^ (l/2)E^(x). 

The generalization of the above expression to the electrodynamic case will be obtained 
simply by replacing in the integrand the electrostatic energy density (1/2)E2 with the full 
electromagnetic energy density = (1/2) (E^ + B^). The same procedure can then be 
applied also with the momentum density T^^ . In this way, for a given inertial reference 
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frame we are led to consider as the electromagnetic energy-momentum contained at the 
time t in a volume V the quantity 

Pi^,.Jt)^'JZ [ d'^T'^{^,t), (2.6) 

where 

pt^u ^ p^.\pu^ _ ^gt^upvxp^^ (2.7) 

is the familiar electromagnetic energy-momentum tensor (we put c = 1 and use the metric 
g^^ = g'^'^ = g^^ = +1, g^^ = —1). Of course, in order for the definition ( p.6| ) to be 
meaningful, the integral on the right-hand side must have a behavior of the form (|2.5| ). It is 
possible to prove that this requirement is actually satisfied as a consequence of the Maxwell 
equations for the electromagnetic field in the presence of point charges. In fact, the field near 
a charged particle can be decomposed as F^^^ = F^^t + F^eu where F/^t is the retarded field 
generated by the particle and F^^t satisfies the free-field equations in a neighborhood of the 
particle. Under the reasonable assumption that F^^t be non-singular in this neighborhood, 
the divergent part of the field can be expressed using the asymptotic expansion for F!^^^ and 
Fret which is provided in App. A. This can be written as 

KL/reti^ + ^) = {Fn-2{X) + {Fn-l{x) + (i^aTvAct)o(x) + 0(x), (2.8) 

with 

(F)„2(a:) = -^«Ax, (2.9) 



— - ) iii ■ x)u Ax — (u ■ x)u Ax u A u 

2r / 2 



(2.10) 



while (Fadv/ret)^, givcu by Eqs. (|A12| ) and ( |A13| ), is a term that is finite for x — 0. In the 
above formulae, z is the space-time coordinate of the particle at a given proper time r . We 
have also introduced the four-vectors u = dz/dr, u = du/dr, the scalar y = ^x'^ + {x ■ u^, 
and we have used the notations u- x = u^x^ , x"^ = x- x = x^x'^ , {u Ax)^'^ = u^x" — u'^x^ and 
similar. One can see that, since ?/ is a homogeneous function of x of degree 1, the term (F)_2 
is homogeneous of degree —2, and moreover it is odd under the inversion x —>■ —x (leaving 
u and it unchanged). Consequently, when the function, as in Eq. (p.6|) , is restricted to the 
three-dimensional space with x° = 0, this term is odd under the space reflection x — x. 
Similarly, (F)_i is homogeneous of degree —1 and is even. It follows that T^'^ has a leading 
term of order which, when substituted to G in an integral of the form ( p.5| ), gives rise 
to the term a'/r on the right-hand side. Besides, the term of T"'^ of order x~^ is odd, so it 
does not generate any logarithmic term in r. One can therefore conclude that the integral 
has indeed the required behavior. 

It can be observed, using the same arguments, that the value of the 7^-integral on the 
right-hand side of Eq. (|2.6| ) is unchanged if the cut-off regions near the point charges, instead 
of being spherical as in Eq. (|2.5|), are chosen with a different shape, provided that such a 
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shape be invariant under space inversion with respect to the particle position. Owing to 
the form of the expansion (|2.8|) -( 2?T0|) , it turns out that mathematical manipulations are 



simpler when surfaces of the form y = constant are considered. Therefore we shall make use 
more often of the equivalent expression 

Pi^^,Jt)=n]imP'y^,Jt,r), (2.11) 

with 

_ r ^ 

Pv,eJt,r)^ / ci3xT°'^(x,t)n%n-r) (2.12) 

•^^ n=l 

and 

yn = \J{.X- ZnY + \Un ' (x - Z^f . (2.13) 

Using = = t, z„ = z„(t), ^^n(t) = dzn/dt, 7„ = (1 — v^) we can write also 



|/„ = V(x - z„)2 + 72 [v„ ■ (x - z„)]^ (2.14) 

from which it can be seen that the surfaces = r are just spheres relativistically contracted 
by a factor 7 in the direction of movement. As shown in App. B, the term of order 1/r in 
Pygjn(t,r) can be easily calculated from the expansion (|2.8|) . This gives 



Pv,eJt, r) = E + Pv,eJt) + 0(r), (2.15) 



Pv,eJt,r) = ^ E ?n {l< - 1%') + PlJt) + 0{r), (2.16) 

where the sums on the right-hand sides extend to all the particles contained in the volume 
V. One can observe from these expressions that the divergent terms of the electromagnetic 
energy-momentum do not transform as a four-vector under Lorentz transformations. In 
particular, for a single particle with f -C 1, the 1/r terms alone would give rise to the well- 

known result Pyem = (4/3)-Pyem^ • However, these divergent terms depend only on the 
particle velocities and contain no information about the interaction of the particles with each 
other or with incident external fields. This suggests that these terms may be dropped as it 
was done for the constant divergent term in Eq. ( p.3|) . In so doing, we are led directly to the 
definition ( p.6|) . In the following subsection we are going to provide a rigorous justification 
for this definition. 



B. The conservation of the renormalized energy-momentum 
and the Lorentz-Dirac equation 

If we denote with Py{t) the total energy-momentum contained the volume V at the time 
t, in a local theory we must expect that 
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^P^{t) = -I dan,r^{^,t), (2.17) 

Clt J dV 

where n is the outward unit vector normal to the surface element da. Clearly the right- 
hand side of the above equation represents the total flow of energy-momentum through the 
boundary dV of the volume V. We are now going to show how the condition expressed by 
Eq. ( p.l7|) , combined with the definition ( |2.6| ) for the electromagnetic energy- momentum, 
univocally determines the dynamics of the charged particles in a relativistically covariant 
way. However, in order to formulate the equations of motion, we first need to introduce the 
fields F'^" (zn), which are defined only on the world-line of the particles according to the 
relation 

F'^^ (zn) = n lim / ^F^^ (z„ + X, t) , (2.18) 

with t = z^. The integral in the solid angle on the right-hand side represents the average 
of the field on a spherical surface of radius r centered at z„(t). Again, one can deduce from 
Eq. ( |2.8| ) and from symmetry considerations that the r dependence of the integral is indeed 
of the form required for the 7^-limit to be properly defined. 

Theorem 1 // one defines 

P^t) ^ Pi^,eJ^) + E rUnKit), (2.19) 

with Py^e^it) given by Eq. ^2. (\ ), then the condition expressed by Eq. ^2. 1 ?| j is satisfied for 
arbitrary V if and only if the particles obey the equation of motion 

du^ 

m„7n^ = qnUn,uF^'\Zn) , (2.20) 



with F^^^ {zn) given by Eq. ( ^.iSj ). 

The total energy-momentum of a system with charged particles is obtained from Eq. 
( p.l9| ) by taking as the volume V the whole space. This gives 

N 

P'it) = PUt) + E ^n<it) , (2.21) 
n=l 

with 

PLit)^nj d^^T'^{-K,t). (2.22) 

From Eq. ( p.l7| ) and from the usual assumption that the field vanishes sufficiently fast at 
infinity, it follows that P^{t) is a constant of the motion, i.e. satisfies 

dP^' 

^ = 0, (2,23) 

Furthermore, one can prove that it has the right transformation property required by special 
relativity. 
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Theorem 2 The total energy-momentum P'^ given by Eq. ^2. 21\ ) transforms as a four-vector 
under Lorentz transformations. 



The proofs of Theorems 1 and 2 are presented in Appendixes C and D respectively. It is 
important to observe that if we have a single charged particle in uniform rectilinear motion, 
and if the field at every point of space reduces to that produced by the particle (i.e. it is equal 
to the retarded or advanced fields, which are in this case coincident), then the expression 
-Pem(^) r), defined by Eq. ( p.l2| ) with V extending to the whole space, is exactly equal to the 
1/r term alone. From the definition of 7^-integral it follows that P^^it) vanishes, so that 
the total energy- momentum is given by P^ = mu^. An analogous result holds for a set of 
uniformly moving particles at very large distances from each other. These considerations 
imply that the parameters m„ appearing in Eqs. (|2.19| )- (|2.2CI|) are to be identified with the 
total physical masses of the particles. 



Considering the equation of motion ( p.20| ), we see that the right-hand side has the fa- 



miliar structure of the Lorentz force, but with the new finite field F^^ instead of the total 
electromagnetic field F^'^ , which of course diverges on the trajectory of the particles. Eq. 
( p.20| ) has the correct mathematical form that an equation of motion in a theory of inter- 
acting particles and fields is supposed to have. It expresses in fact the second-order time 
derivative of the particle's position as a function of the position, its first-order derivative 
and the instantaneous field configuration, since F^^{zn) is a function of the field F^'^(x, t) 
at fixed t in a neighborhood of Zn{t). A complete description of the temporal evolution of 
the system is obtained by coupling with Eq. ( |2.2U| ) the Maxwell equations 

^?^ + VAE(x,t)=0 (2.24) 
- V A B(x, t) = -j(x, t) = - E QnMt)S (x - z„(t)) , (2.25) 

n=l 

with the supplementary equations 

N 

V ■ E(x, t) = p(x, t) = Y^ qj (x - z„(t)) (2.26) 

n=l 

V-B(x,t) = (2.27) 

to be regarded as usual as constraints on the initial conditions. 

In order to clarify the meaning of Eq. ( p.2(]| ), let us define for a given charged particle 
the field 

F'' ^F'^'-l iKL + FZ) , (2.28) 

which remains finite on the world-line of the particle. Using Eqs. ( |2.8| ) and ( |A12| ) we can 
then write 

F'^'^^z + x) = (F'^")_2 {x) + (F^^").! (x) + (F'^")o (x) + F^\z) + 0{x) , (2.29) 
where {F^^)q{x), given by Eq. (|A13D, is a term that changes sign under the substitution 



X —X. Therefore, applying the definition ( p.l8| ) and taking into account the spatial 
antisymmetry of {F)_2 and (-F)o, we obtain 
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F'"'(z) = lim 



da 



X 



OJ|x|=r An 



{F^'')^{x) + F^\z) =F^''{z) (2.30) 



(a proof of Eq. (|2.30| ) for the system of reference comoving with the particle was already 
given in Ref. 0). Since the definition ( |2.28| ) is Lorentz covariant, Eq. ( p.30| ) shows that 
F^'^(z) actually transforms like a tensor. Substituting F^" to F^'^ on the right-hand side of 
Eq. ( |2.20D , one can recognize in the latter a well-known form of the relativistic Lorentz-Dirac 



equation. We note also that Eqs. (|2.8| ) and ([A12|) are in agreement with the relation first 
proved by Dirac in []T| 

I - i^al(^)] - iu'u' - n^nn for x 



from which it follows that, after decomposing the field as F^'-' = F^^t + -^/et > the Eq. (p.20|) 
can be rewritten as 

mii^' = qu.FZ + l^iuf^ - iiuU^u^, (2.31) 
3 47r 

where use has been made of the relations u"^ = —1, = d{u'^)/dT = 2u ■ u, = d{u ■u)/dT = 



u 



2 



+ il ■ u. In the familiar form (|2.31|) the Lorentz-Dirac equation appears explicitly as 



a third- order differential equation allowing to calculate the trajectory of a point charge 
interacting with an assigned "external" field Fg^t (given in general by an incident radiation 
field plus the sum of the retarded fields generated by other charges). On the contrary, in our 
formulation of the theory Eq. (|2.20| ) is a second-order equation which alone cannot determine 
the particle's trajectory, since it is intrinsically coupled with the Eqs. ( ^.24[ ) and ( ^.251 ) for 



the total field F^'^ . The three equations need be solved together in order to determine 
the temporal evolution of the system of interacting particles and fields. One may observe 
incidentally that these equations are manifestly invariant under time reversal, which is not 
the case for the Lorentz-Dirac equation written in the form (|2.31|) . 

The solution of the system of coupled equations ( |2.20| ) and ( |2.24| )-( p.25| ) depends on a set 
of initial conditions that is entirely specified by the values, at a given time t, of z„(f), v„(f) 
and F^'^[yi,t) for every x. These initial conditions are subject to some relevant constraints, 
in addition to those already expressed by Eqs. ( p.26| ) and (|2.27|) . The reasonable requirement 



that F^'^{x) be non-singular at Zn implies that the divergent part of F'^'^(x,t) must have 
the form given by Eqs. ( p.8|) -( prT0D . Therefore the term of order |x — z„(f)|~^ is completely 
determined by v„(f). Furthermore, the term of order |x — Zn{t)\~^ is related to v„(f) and to 
the zeroth order term F'^'^ (z„(f)), since the vector u„ appearing in Eq. ( p.lQ ) must coincide 
with that resulting from Eq. ( p.20|) . 

Another subtler constraint on the initial conditions is related to the well-known problem 
of the "run-away" solutions of the Lorentz-Dirac equation [|l]-|^ . We remark incidentally that 
such solutions do not violate the conservation law given by Eq. (|2.23|) . In fact the unlimited 
increase towards -|-oo of the particle's kinetic energy is compensated by a corresponding 
decrease towards — oo of the electromagnetic energy Pg*]^ (which is not a definite positive 
quantity, according to our definition), and analogous statements are valid for the three spatial 
components of the momentum. These solutions, which are therefore completely legitimate 
from a mathematical standpoint, are generally dismissed for phenomenological reasons, since 
run-away electrons have never been observed. 
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We remind that, in the usual treatment of the Lorentz-Dirac equation as a third-order 
differential equation, the initial particle acceleration has to be specified among the initial 
conditions, together with the velocity, the position and the incident field. One may then take 
advantage of this additional free parameter in order to get rid of the run-away solutions. For 
instance, considering a scattering problem between a particle and a wave-packet initially very 
far apart, it is possible in general to obtain a finite result for the asymptotic momentum of 
the particle in the final state by introducing a suitable "preacceleration" , i.e. an acceleration 
of the particle that decreases as exp(t/7ro) for t — > — oo, with tq = q'^/Gnmc^ (~ 0.62 x 10~^^ 
sec in the case of an electron), and that becomes significant starting from a time of order tq 
before the actual interaction between particle and wave-packet takes place. 

In our formulation of the Cauchy problem for the equations (|2.20|) and ( p.24|) -( p.25D , 
since the particle equation of motion is now of second order, the initial acceleration cannot 
be independently assigned, but is determined by the initial value of the total electromag- 
netic field according to Eq. ( 2.20 ). Therefore the requirement that the asymptotic particle 
momenta be finite both for t — oo and for t +00 is equivalent to a constraint on 
the initial values of the field. The fact that, in the traditional "third-order" treatment, the 
preacceleration is exponentially vanishing for times t very far in the past with respect to 
the actual interaction, within the new "second-order" picture translates into the fact that 
an equivalently small "adjustment" of an arbitrarily assigned field configuration at the same 
times t is generally required in order to satisfy the asymptotic conditions. 

It is finally worth mentioning that the run-away behavior cannot be revealed by a per- 
turbative expansion in the coupling constant. Therefore the existence of run- away solutions 
cannot be claimed to be a formal drawback of the classical theory of point charges as 
compared to perturbative quantum electrodynamics. It has also been conjectured that a 
particular kind of run-away solutions may exist in relativistic spin 1/2 QED 



III. THE ACTION PRINCIPLE 
A. The renormalized lagrangian 



A fundamental property of the system of differential equations constituted by Eqs. (p. 20 ) 



and (p.24|) -( p.27|) is that it can be derived from an action principle. A lagrangian formalism 



can in fact be introduced in a natural way by taking as configuration variables of the system 
the position z„(t) of the particles and the four- vector electromagnetic potential A^(x, t) 
satisfying the relation 

Given an inertial system of coordinates, we shall define as the total lagrangian of the system 
at the time t the functional 

L = i^em + -^part + -^int ; (3-2) 

with 

Lem = -\n I c/=^xF,,,(x,t)F^'^(x,t), (3.3) 
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^ m 

Vrt = - E — ' (3-4) 

n=l 

i^int = E -"n^'^ (^n) ' (3-5) 
n=l 7" 

in Eq. ( p.5| ) is defined in perfect analogy witli Eq. ( p.l8| ): 

A" = 7^ lim / (z„ + X, t) . (3.6) 

One can see tlierefore tliat L is indeed a functional of the configuration variables and of 
their first time derivatives, taken at a given time t. In general this functional is defined for 
every field configuration having a behavior near the singular points such that the 7^-integral 
in Eq. (|3.3| ) and the 7^-limit in Eq. (|3.6| ) are properly defined. 

As usual, Eqs. ( |2.24| ) and ( p.27| ) are automatically satisfied as a result of Eq. ( |3.1| ), so 
that the equations of motion for the variables z„ and A^(x) reduce to Eqs. ( 2.2CI| ), ( p.25|) and 
( [^.2(j| ). We are now going to show that these equations follow from Hamilton's variational 
principle. 

Theorem 3 // one defines the action I as 

rt2 



L dt, 

ti 

with L given by Eqs. 3^ )-l ^~3i ), then the equations of motion of the system in the time 
interval ti < t < t2 are equivalent to the condition that 61 = for every infinitesimal 
variation 6A^^{x) and 6zn{t) vanishing at t = ti and t = t2, and such that the varied 
trajectories {z„ + 5z„,yl^ + 6A^} remain in the definition domain of I. 

In order to prove this theorem, we shall introduce the non-renormalized lagrangian 

L(r) = Lem(r) + Lpart + Ant('^), (3.7) 



with 



and 



- 1 r ^ 

LeUr) = -T / d='xF^,(x,t)F^^(x,t) [] ^ (|x - z„| - r) (3.8) 

^ n=l 



LUr)-Z;f<lJ-^M^.^^,t). (3.9) 

Since L = TZlimr-,o L{r), we can write 

L(r) = (L)_ir-^ + L + 0(r) , (3.10) 

where the coefficient {L)-i is some other functional of the configuration variables and of 
their first time derivatives. It follows from Eq. (|3.1CI|) that 

6L{r) = 6{L)_ir-^ + 6L + 0{r) . (3.11) 

Therefore the variations of the renormalized lagrangian L will be evaluated by calculating 
the corresponding variations of L, and then taking their 7^-limit for r — 0. 
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B. Derivation of the field equations from the action principle 



Let us fix at a first stage 5z„ = and consider only variations of tlie electromagnetic 
potential, which will be represented by arbitrary infinitesimal test functions SA^^lx) satisfying 
the boundary conditions. We want to show that, in order to have identically 61 = 0, the 
fields F'^'^ must define distributions that satisfy the differential equations ( |2.25|) and ( p. 261) . 
For 6z„ = we have 



5L{r) 



r ^ a 

/ d'^F^^'^d^dA, + J2 —u'JA^iz^) + O 

n=l '^n 



(3.12) 



where we have used the fact that, assuming F^*^ to have at most a singularity of order r ^ 
in Zn-, the integral on the right-hand side is convergent, and moreover 



dVLy 

|x|=r An 



■5A^ (z„ + X, t) = 5A^{zn) + 0{r^). 



Since on the right-hand side of Eq. ( p.l2| ) there is no term of order r ^, a comparison with 
Eq. ( p.ll|) shows that 5{L)^i = and 



61 = r dt6L = r dt \- f d^xF'^'dJA^ + V —u'^^SA^ 

Jti Jti [J n^lln 



(3.13) 



Since SA^, are test functions vanishing for t = ti and t = t2, Eq. (|3.13|) can be rewritten as 



SI = j^J dt / d^x 



N 



n=l 'fn 



6A„ 



where the partial derivatives of F^'^ are defined according to the distribution theory. There- 
fore the condition 61 = for arbitrary SA^y implies 



N 



qn 



n=l /" 



Z,„ 



(3.14) 



which is equivalent to Eqs. ( |2.25| ) and (|2.26| ). Of course, assuming that the fields F^^ are 
represented by ordinary differentiable functions of x everywhere outside the world-lines of 
the particles, Eq. ( |3.14| ) implies in particular that in the same region such functions satisfy 
the differential equation d^F^^^ = in the usual sense. 

It is well-known that the general solution of Eq. ( p.l4| ), when the particle trajectories are 
assigned, is given by the sum of the Lienard-Wiechert fields and a solution of the free-field 
equations. Using this fact, we are now going to show that the lagrangian L is indeed correctly 
defined on the solutions of the field equations and on the neighboring trajectories obtained 
via non-singular variations SAj^ of the potential. This fact ensures the self-consistency of the 
whole procedure. First of all, since asymptotic conditions of the type (p^ ) must be satisfied 
by the fields near the particles, we can repeat the same considerations we made about the 
7?.- integral in Eq. ( |2.6| ), showing that also the right-hand side of Eq. ( p.3| ) is correctly defined. 
Moreover, using the expansion for the advanced and retarded vector potentials derived in 
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App. A, one can verify that also the potential A„ appearing in Eq. (|3.5| ) is well defined 
according to Eq. ([3.61). In particular, defining for a given particle 



— ( Af^ -1-4'^ 

^ — ^ 2 adv + ^re1 



one gets from Eqs. (^) and (|A8D 

A^^iz + x) = (A'^)^i (x) + (A'^)o (x) + A^{z) + 0(x), 
and from the antisymmetry in x of {A'^)q (x) it follows that 



A^iz) = lim 



Oj|x|=r 4:71 



{A^o (x) + A'{z 



A'(z) 



(3.15) 



(3.16) 



(3.17) 



in close analogy with Eqs. ( p.28| )-( p.30[ ). We remark also that Eqs. ( p.l5| ) and ( p.l7[ ) ensure 
that A^{z) behaves as a four- vector under Lorentz transformations. 



C. Derivation of the particle equation of motion from the action principle 

When considering also variations 5z„(t) 7^ 0, one has to pay attention to the fact that the 
quantity A,y (z„) is defined only on the world-line of a particle, so that its partial derivatives 
with respect to z„ are not defined in the usual sense. This is equivalent to the obvious 
statement that, whenever 5z„(t) 7^ and 6A'^{x) = 0, the varied trajectories no longer 
belong to the definition domain of the functional L. In order to remain within this domain, 
each variation 6zn(t) 7^ must be accompanied by a simultaneous variation 6A'^ of the 
potentials, characterized by an appropriate singular behavior near the particles. For this 
reason, it is convenient to proceed in two steps. First, given the boundary conditions at 
t = ti and t = t2, for each possible trajectory z^ of the particles one determines the potential 
A'^ (x) that minimizes the action / with respect to the variations of the potential alone: as 
we have seen in the previous subsection, ^^^(a;) is just the potential that solves the Maxwell 
equations for those particular particle trajectories. Then one can define the functional I' 
of the trajectory of the particles alone, as 

l'[Zn]=l[zn,Al{x)] . (3.18) 

It is now possible to show that the condition 61' = 0, for arbitrary variations 5z„(t) satisfying 
Szniti) = 6zn(t2) = 0, is equivalent to the fact that z„(t), together with A^^(x), are solutions 
of Eq. (|2.2CI|) . Of course, the minima of / must coincide with the minima of /': therefore the 



trajectory of particles and fields |z„,A^^(x)| satisfying Eq. ( p.20|) is indeed the trajectory 
that minimizes / within its whole definition domain, subject to the assigned boundary 
conditions. 

Keeping these ideas in mind, we shall still write for simplicity A'^ instead of A'^^ , although 
we shall always implicitly assume that the field satisfies Eq. (|3.14 ). We shall therefore apply 



the previous formulae describing the asymptotic behavior of the field near the point charge. 
We shall make use also of an expansion for the derivatives of the potentials that, similarly 
to Eq. (|2.29|) , follows from the results of App. A: 
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d^A,{z + x) = {x) + {x) + {d^A,)^ {x) + d^A,{z) + 0(x), (3.19) 



where {d^A^). {x) is a homogeneous function of degree j and even (resp. odd) parity for j 
odd (resp. even), while 



duAJz) = 7?.hm 



d^Ay (z + X, t) . 



(3.20) 



Note that the above equation is to be considered as a definition since, as we aheady observed, 
the partial derivatives of A{z) do not exist in the usual sense. However one can easily show 
that, in close analogy with Eqs. ( p.30|) and ( |3.17| ), 



d^A^{z) = di,A^{z), 



(3.21) 



where the right-hand side represents the usual partial derivative of the function A{x) defined 



by Eq. ([3151) . 

With a rather tedious but straightforward calculation, which is reported in App. E, one 
obtains 

d 



N 



+ ^ \(inAi{zn)5v\ + g„7„ ^v!'^diA^{zn)6z\ + miUn,i5v\ 
=1 

- N _ 

^ (m„M„,i + qnAiiZn)) Szi,-n d'l^F^^A, 

.n=l 



n=l 

dt 



N 



-E 

n=l 



d 

dt 



(mnUn,i + qnAi{zn)) - g„7„ 



Szt 



(3.22) 



Since = and 6Ai, = for t = ti and t = t2, the boundary terms do not contribute to 
the variation of the action, and we obtain 



SI = - f dtJ2 4: {^nUn,i + qnAi{Zn) \ - Qnln ^KdiA^{Zn) 
n=l L'^'' 

Since it follows from Eqs. (|3.17|) and (p.21|) that 



Szt. 



(3.23) 



the condition 61 = for arbitrary 6zl^ implies 

du 



'"n,i 



dt 



qn 



{diAy{zn) - dyAi{zn)) < = qnFi^{zn)u: 



V 

n 1 



(3.24) 



that coincides with Eq. ( p.20[ ). We have therefore proved that the dynamics of the system 
is entirely determined by Hamilton's variational principle applied to the lagrangian defined 
by Eqs. (P7^)-(P75|), as stated by Theorem 3. Of course, it is also possible to deduce directly 
the conservation laws of energy and momentum from the translational invariance in space 
and time of L. This is illustrated in detail in App. F. 
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D. A comparison with the idea of mass renormalization 



The renormalization of our lagrangian was achieved simply through the introduction of 
a 7^-integral for the electromagnetic part in Eq. ( |3.3|) and a 7?.-limit of the electromagnetic 
potential for the interaction part in Eq. ( pj.SI ). Clearly this procedure has in principle nothing 
to do with the renormalization of the particle's mass, which is typical of other approaches 
to the argument, and in particular of standard renormalization techniques of quantum field 
theories. Therefore no infinitely negative "bare mass" has to be introduced in our theory, 
and the difficulty of simultaneously compensating with such a mass term the non-covariant 
divergent parts in Eqs. ( p.l5[ )-( p.l6[ ) is automatically avoided. It is however interesting to 



-ir 



note that the divergent part of Eq. (| 

1 

4 
1 

47 



^) can be easily calculated using Eq. 
2g2 



-— y'dy da 



(47r)V 



(we consider here for simplicity a single charged particle), so that we have 

1 



lim 



--J d'^F'^^F^Jiy 



q- 



(3.25) 



Similarly we have from Eq. 



so that 



-u'^n hm 

'J r— >0 



dfl^ 



An 



d^y (A^) 



dQ 



lim 



7 



dQ^ 
Air 



r 



■AJz + x) + 



Therefore, combining Eqs. 



lim 



-- / d^^F^^'F, 



0{y 



and ( p. 261 ), we can write 



+ -u 

7 



y=r 



dQy 



.AJz + ^)-^ 



mo 

7 



(3.26) 



(3.27) 



with 



mo 



m 



Clearly, one could interpret Eq. ( |3.27D as the lagrangian of a particle having indeed a 
divergent negative bare mass mo and carrying a total electric charge q that, for an observer 
with the same instantaneous velocity of the particle, looks uniformly distributed over an 
infinitesimal spherical surface of radius r. We must however underline that Eq. ( 3.27|) is 
verified only when the fields are already solutions of Maxwell equations. In general, taking 
Eq. (|3.27| ) as a definition would not be equivalent to the definition of L that we have given 
in Eqs. (|3.2| )- (|3.5| ). This becomes particularly evident when introducing the hamiltonian 
formalism. In fact, since each of the three terms on the right-hand side of Eq. ( |3.27| ) is 
divergent, it is immediate to see that also the conjugate momenta of the particles Pi = dL/dv^ 
are divergent in the r —>■ limit. On the contrary, the lagrangian { \i.'2\) is the sum of three 
individually finite terms. As we shall show in the next section, from this fact it follows that 
all the momenta are finite, and therefore an exact hamiltonian can be defined for the theory. 



16 



IV. EULER-LAGRANGE AND HAMILTON EQUATIONS 



It is useful to define the partial derivatives of the lagrangian L with respect to its ar- 
guments, in such a way that the equations of motion can be formally derived as the usual 
Euler- Lagrange equations associated with L. We have seen in the previous section that, in 
order to remain within the original definition domain (that we shall call Dq) of L, each vari- 
ation 6zn must be accompanied by a variation 6 A of the electromagnetic potential with an 
appropriate singular behavior near the particles. Therefore, defining a partial derivative of 
the type dL/dzl^ for constant potential A^{x) requires introducing an extension of L outside 
-Do- We have to note first of all that, inside Dq, the particles' positions Zn are univocally 
determined by the field variables A^^{x), since the potential is assumed to be singular just in 
Zn- In other words, as soon as y4'^(x, t) is assigned for every x, also its singular points, that 
will be now called w„(t) to distinguish them in general from the particles' positions z„(t), 
are of course determined, and inside Dq it is always assumed that z„ = w„. It is clear now 
that, in order to treat z„ and A^^{x) as really independent variables, we have to define L 
also for z,„ 7^ w„. Having this in mind, the results of the previous section indicate how to 
extend in an obvious way the definition of L. We shall put 



dzi 



0, 



(4.1) 



which is consistent with considering as the purely electromagnetic part of the lagrangian, 
and 



dL 



o • — (ln~diAi^(Zn) , 



(4.2) 



which is equivalent to extending in a neighborhood of w„ the function A^{zn) appearing in 
Eq. (|3.5|) according to 



A^(z„, t) = A^(w„, t) + diA^iwn, t)(4 -wi) + 



. Zr, - W.„ 



On the right-hand side of the above equation the quantities A^{wn,t) and diAj,{wn,t) are 
defined according to Eqs. ( p.6|) and ( |3.2CI| ) respectively, whereas terms of order (z„ — w„)^ 
can be neglected for our purposes, since they clearly do not affect the equations of motion. 
Besides, the derivatives with respect to v^^ = dzl^/dt can be determined unambiguously as 



dvi 



0, 



qAi{zn) 



We have then, as a result of Eq. ( p.24|) 



Qn t)iA^[Zn) = — TTlnUnA + QnAiyZn) 

at 



dzl 



In 



d dL 

dt dvl 



Similarly, using Eq. (|3.14| ), we have 
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5L 



N 



5A„{x) 7„ dtddQA^[x) 



where the functional derivatives of L with respect to A^{x) and OqA^^x) are unambiguously 
defined as distributions in the usual way. Looking at the first and last members of the two 
above equations, the formal structure of Euler-Lagrange equations is directly recognizable. 

The hamiltonian functional can be introduced via a Legendre transformation on L, as in 
ordinary analytical mechanics. As usual, after imposing the Lorentz gauge condition 

d^A^{x) = 0, 

we shall add a term —(1/2) (S^A^')^ to the electromagnetic lagrangian density in order to 
introduce a conjugate momentum also for the temporal component A^{x) of the potential. 
We shall therefore replace Eq. (|3.3|) with 



As a result, the Euler-Lagrange equation for A^ becomes 



N 



d.d'^A^'ix) 



n=l 



-fix) 



Then, after introducing the canonical momenta 

dL 



Pn,i 



dvi 



6L 



TlAX = 



noix = 



6dnAHx) 



6L 



d,A^{x), 



we can define the hamiltonian H of the system as 

N 

H = Y^ Pn,K + 7^ / rf^x vr^^o - L 

n=l 

= njd'^ I^^FijF'^ + ]^'K^'K^ - 7^'diA^ + 7r°9iA' 

N 

+ E 

n=l 



+ ^j {pn,i - qnAi{zn)) {pi - qnA'iZn)) + ml 



Hamilton's equations are: 



dzl 



dH 



pI - q„A'{z„ 



dt dpn,i 



\j {Pn,k - QnAkiZn)) (p^ " qnA''{Zn)) + ml 



(4.3) 

(4.4) 
(4.5) 
(4.6) 



(4.7) 



(4.8) 



doA\x) 



SH 

6TTi{x) 



n'{x)+d'Ao{x), 



(4.9) 
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doA'^ix 



5H 



n\x)-d,A\x), 



dPn 

dt 



do7r'{x) 



dH 



5no{x 

qnd'A^{Zn) {pn,j " QnAji 



6H 



5Ai{x 



\j {Pn,k - QnAkiZn)) {p^ - qnA^{Zn)) + m. 

qn{Pn- qnA'iZn)) 6{yi~ Zn 



qnd'A^iz^), 



(4.10) 
(4.11) 



n=l 



{Pn,k - QnAkiZn)) (p^ - g„A'=(z„)) + 



9o7r°(x) 



5H 



N 



5An{x) 



di7i'{x) + g„5(x - z„). 



(4.12) 



(4.13) 



^ol-y „=i 

As usual, Eqs. (^^ ) -( ^?ToD are just the inversion of Eqs. ( [4.4|) -( ^^ , whereas Eq. ( f4.11|) is 
equivalent to Eq. and Eqs. (^1^) - (^A3|) are equivalent to Eq. ( ^ . 

The theory can be expressed by means of fundamental Poisson brackets between canon- 
ical variables having the familiar form 

{A'^(x),7r'^(y)} =^'^''<5(x-y), (4.15) 
all the others being equal to zero. After introducing the total momentum 

N 



n=l 



and putting P° = H, it is straightforward to verify that from Eqs. 

apo dzi 



(4.16) 
-(^J^ it follows 



{Pn,P'} 

{A>^ix),P'^} 
{n>^{x),P''} 



dpn,i 
dpn,i 

dP^ _ 

dZji^i 

SP" 
5tt^{x) 
SP" 
SA^{x) 



dt 



dpi 
dt 







-d'A^ix) 
-- -d'^Ti^'ix) 



where the members of the last two equations are to be interpreted as distributions. The 
above equations show that, as expected, P'^ is the generator of the space-time translations 
for the system. Of course, on the solutions of the equations of motion it is equal to the total 
energy-momentum introduced in Sec. II. This is explicitly shown in App. F. 
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V. THE RENORMALIZATION PROCEDURE IN THE WAVE- VECTOR SPACE 



A. The energy-momentum in terms of the Fourier transform 
of the electromagnetic field 

Let us operate the Fourier decomposition of the electromagnetic potential according to 
the formulae 



1 



rf^k r 



a^(k,t)e'''''' + a*(k,t)e 



-ikx 



2(27r)= 



a^(k,t)e^^-''-a;(k,t)e- 



ikx 



(5.1a) 
(5.1b) 



with = |k|. As a consequence of Eqs. ( ^.15| ) and ( p.l[ ), the Poisson brackets between the 

{a^(k), a.(k')} = {a;(k), ^(k')} = (5.2a) 



variables and a* are 



{a^(k), <(k')} = ~2{2'Kfig^,kH{\i - k'). 



(5.2b) 



It is well-known that for a free field one can express the electromagnetic energy- 
momentum as 



2(27r)3 



Fa'^(k)a*(k). 



It is not surprising that, when point-like charged particles are present, also the above integral 
becomes divergent. The same is true for the integrals on the right-hand side of the Fourier 
expansions ( p.l| ) when x takes the value of the position z„(t) of one of the particles. This 
can be verified using the following asymptotic expansion of a^{k,t) for |k| oo, which is 
derived in App. G: 



N 

a'^(k, t) = Y. [«)_! (k, t) + «)o (k, t)] e-^-- + O (|kr^) 



n=l 



with 



(an)_i 



Ur. 



-iq, 



k ■ Uji 

{k ■ Un)Un - {k 
{k ■ Unf 



(5.3) 

(5.4) 
(5.5) 



Clearly, Eq. ( ^.3|) is the direct counterpart of the expansions (|3.16| ) and ( p.l9|) of A'^z + a;) 
and df^A'^{z + x) for |x| 0. We have in particular 



(A^n), (^) 
(9.A^)_ (x) 



2(27r)3 7 fco 
i r d^k 



<),(k,t)e^'^-'' + «);(k,t)e 



— ikx 



2(27r)3 7 k^ " 



t:)^.(k,t)e^'^--«);(k,t)e- 



-ikx 



(5.6) 
(5.7) 
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It can be noticed that (cin)j is a homogeneous function of k of degree —j — 2. Besides, it is 
real for j = —1 and purely imaginary for j = 0. 

We can however redefine the divergent integrals in k space by means of a renormalization 
procedure equivalent to that introduced in Sec. II for the corresponding expressions in spatial 
coordinates. 

Definition 2 Given a function /(A) having for A — > cxd a behavior of the form 

/(A) = cA + ci + 0(A-i), 
we define its TZ-limit for A ^ oo as 

7^ lim /(A) = d. 

A^oo 

Correspondingly, given a function g{k) having a behavior such that 

I rf^k ^(k) = c'A + d! + Oik-^) for A ^ oo, 

J|k|<A 

we define its TZ-integral over the k space as 

n f d^kg(k) = 7^ lim / d^kg(k) = d! . (5.8) 

J A^ooJ|k|<A 



Using this mathematical tool one can express in terms of the variables and a* the 
canonical energy-momentum four-vector for the system of interacting particles and field, as 
stated by the following theorem. 



Theorem 4 The four-vectors and A^^(zn), given respectively by Eqs. ( \4-V , o,nd 



(^Tdj), can be equivalently expressed by the equations 

d^k ^ 



and 



P 



1 r d-^k 

^^7^ j — A;V(k,t)a:(k,t) + X:^^ 

P' = ,^^n j d'ka^{k,t)al{k,t) 

N 

+ E 



n=l 



^/i(^n) 



qnA^{Zn) + \j {pn,i - qnAi{Zn)) (K " gn-4*(2;„)] + 



2(27r) 



^7^ 



d^k r 



A;0 



a^(k,t)e''^"" +a*(k,t)e 



(5.9) 



(5.10) 



(5.11) 



The proof of Theorem 4 is given in App. H. By making use of Eqs. (5^), ( 5.10|) and 



( p.ll| ), one readily obtains the equation of motion of the 



21 



d 
It 



(k,t) = {a'^(k,t),P°} = -zA;V(k,t) 



2(27r)3A;0„^,^"7n 



?A;V(k,t) + 



2(27r)3A;0 



- jkx 



which is of course in accordance with Eqs. (|4.3| ) and (5J.). We remark also that, with exactly 
the same procedure followed to derive Eq. ( |5.11| ), one can prove that the quantity dyA^{zn) 
defined by Eq. ( |3.20| ) can be expressed as 



2(27r)' 



^7^ 



A;0 



B. Another remark on the renormalization of mass 

We saw in Sec. II that the expression for the electromagnetic energy-momentum in x 
space can be most conveniently handled when the cut-off regions in the integration domain 
are taken as ellipsoids like in Eqs. and ( |2.12| ). Similarly, in some cases the 7?--integrals 



in k space appearing in Eqs. (^^ ) -( ^TToD may be more suitably expressed by choosing a 
different cut-off instead of the sphere of radius A introduced in Eq. (p. 81). A particularly 
interesting situation occurs when a single charged particle is present. In such a case, putting 



with = |k|, it can be easily seen that a'^ transforms as a four- vector, as it can already be 
argued from the form of the right-hand sides of Eqs. (|5.4| ) and i f).5\j . It follows that also the 
expression 

^em(i,A) = y — A;V(k,t)<(k,t)^(A + fc,«^(t)) (5.12) 

transforms as a four-vector, since it can be rewritten as 

KUt,^) = I d^kk^lf{k,t)lf^{k,t)5{e)e{ky[k + k^u\t) 



We observe that the integration domain in Eq. (|5.12D is characterized by the condition 



Ikl - k ■ V < — , 

7 

which corresponds to the volume contained into a rotational ellipsoid with a focus at the 
origin. Repeating the procedure used in App. H to obtain Eqs. ( |H10| ) and ( [Hll| ), it is easy 
to see that 

7^2im irL(t,A) = j — W(k,t)<(k,t) . 

Furthermore, for v(t) = we get from Eq. 
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1 



2(27r)3 

-,2 



q 



2{2n) 



An' 



Thanks to the Lorentz covariance, we can immediately generahze the above result for arbi- 
trary v(t), obtaining 



!(2 



(5.13) 



A similar procedure can be followed for the 7^- integral that expresses Af^{z) according 
to Eq. (|5.11| ). Applying again the covariant cut-off we can define the four- vector 



2(27r)3 
1 

(2^ 



fk 



a^(k, t)e^^-^ + a;(k, t)e-'^-^] (A + k^u^it 

d^k [a^{k, t) + a*^{k, t)] 6{k^)e{ky (A + kxu\t) 

Then, proceeding as in Eqs. (|H17|) -( ^18|) , it is immediate to see that 

A^{z) = 7^ lim A^{z,A). 



(5.14) 



A^CXD 



Moreover, since for v(t) = 



1 



2(27r)3 
9 

(27r) 



A;0 



, (k, t)e*'^-^ + {a,Y_, (k, t)e-*'^-^ 9 A + kxu\t) 



d-'k 



■(A:0)^ 



9 A-k^ 



2n' 



rAM^(t), 



the Lorentz covariance implies that for arbitrary v(t) 



Af,{z,A) 



271- 



Mt^it) + Af,{z) + 0{A-^) for A ^ oo. 



(5.15) 



Putting together Eqs. (|5.9| ), (|5.10| ), ( ^.13| ) and (|5.15|) , we can finally write 



n 



lim 

A— »oo 



2(27r)3 
mQU^'{t) + qA^'{z,A) + 



d^k 
1 



k^'a^{k,t)al{k,t) 
d^k 



2(27r) 



k^ 



fc'^a^(k, t)<(k, t)^ (a + kxu\t) 



with mo = m — (g^/47r^) A. We obtain therefore a result very similar to that discussed in Sec. 
Ill when a proper cut-off is chosen and only fields that already obey Maxwell equations 



are considered, the method of the 7^-integrals leads to expressions that formally correspond 
to the introduction of an infinite renormalization of the particle's mass. The analogy with 
the situation occurring with covariant regularization methods in quantum field theory is 
evident. It will be analyzed more closely in the following section. 
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VI. THE COMPARISON WITH QUANTUM ELECTRODYNAMICS 



In the previous sections of this paper we have provided a complete diverge-free formu- 
lation of the classical electrodynamics of point charges. In the light of this result it is of 
course extremely interesting to reconsider the problem of the relationship between classical 
and quantum theories. It is in fact well known that the difficulties related to ultraviolet 
divergences in quantum electrodynamics have been up to now circumvented only in the 
framework of a perturbation expansion in the coupling constant a = e'^/Airhc ~ 1/137. 
Finite non-perturbative expressions for the equations of motion or the hamiltonian operator 
have not yet been established, and it is even doubtful that they actually exist within the 
boundaries of pure electrodynamics. In connection with this problem, the implications of 
our new way of formulating the classical theory can be directly investigated by making use 
of the correspondence principle. Of course, in the quantum theory one encounters also the 
class of divergences associated with vacuum polarization, that have no direct equivalent in 
the classical theory. We shall not here take them into consideration, and we shall limit 
our analysis to the simplest situation in which the analogy with the classical case can be 
significant, i.e. the electron self-energy to the lowest order in a. We shall assume that the 
charged particles are represented by spin 1/2 fermions (electrons), and therefore described 
as the quanta of a Dirac field tp. We shall follow the scheme of canonical quantization and. 



by analogy with the classical Poisson brackets ( |4.15| ) and ( |4.14| ), we shall postulate the usual 
fundamental equal-time commutation or anticommutation relations between field operators. 

Considering the usual Gell-Mann and Low perturbative expansion of the two-point 
fermionic Green function iS'p{x) = (0|T ip{x)ip{0) |0), the contribution of order a (using 
natural units with c = h = 1) can be expressed as 



+00 



dh / dh (0|T Hi{h)Hj{t2Mx)ilj{0) |0)(o) 



(6.1) 



where the subscript (0) indicates that the matrix element involves the non-interacting vac- 
uum state and field operators, and only connected Feynman diagrams have to be taken into 
account. By analogy with Eq. (|4.7|) or (|5.10|) , the interaction hamiltonian is obtained by 
substituting the field A'^ in place of A^^ into the familiar expression of the non-renormalized 
theory. This gives 



Hiit) 



d^xAf,{x)j'^{x) = -ie / c/^xA^ : ^p-f''^jJ 



(6.2) 



where A^ is given by a 7^-integral of the form ( ^.11| ). The natural way to generalize our 
notion of 7?.-integral when the integrand is an operator in a Hilbert space is to perform the 
limiting procedure in the weak sense, i.e. on the matrix elements of the operator between 
any pair of states. This directly provides the prescription for calculating iSp . One first has 
to replace the field A^(x) in Eq. ( |6.2| ) with the expression 



A^{x,A) 



2(27r)^ 



d'k r 



A;0 



a^(k)e*'=-^ + <(k)e-*'=-^' 6 {A ~ /(k)) 



(6.3) 



where /(k) is an appropriate function of k. Then, after taking the expectation value on the 
non- interacting vacuum state of the time-ordered product indicated in Eq. ( |6.1|) , one has to 
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evaluate a properly defined 7^-limit of the resulting expression for A oo. One may notice 
the direct correspondence between Eq. (|6.3|) and Eq. (|5.14| ) when putting /(k) = —k\U^ 
and considering that, since the field operators in Eq. ( |6.1| ) satisfy the free-field equations, 
one can make the substitution 



a^{k,t) a^(k)e" 



It is easy to see that as a result of these prescriptions the expectation value in Eq. ( p.l\) 
has to be calculated by making use of the following free propagators: 



(0|r[A^(x,AK(0,A)]|0)(o)=2^7^'^ 
{0\TU{x)m] |0)(o) = z 



e'''-''DFik,A) = ig^"'DFix,A) 



e'P-^'SF{p)=zSF{x), 



with 



DFik,A) 



A;2 + A2 - le 



9{A-f{k)) 



and 



Sf{p) 



ip ■ 'J — m 

p2 _|_ ^2 _ ^g- 



As usual, a photon mass term A has been included to remove infrared divergences, but 
it is expected that the limit A — > can ultimately be taken for all physically significant 
quantities. By applying Wick's theorem one finds for the Fourier transform of Eq. (|6.1| ): 

sPip) ^ J d*xSP{x)e-'P-' = Sf{p)^^'\p)Sf{p), 



with 



and 



=7^ lim S^'Hp^A) 

A^oo 



(6.4) 



S(2)(p,A) = d''kDF{k,A)^^SF{p-k)r. 

In the last equation, the integration in dk^ can be performed by closing the integration path 
in the complex plane and applying the residue theorem. This gives 



S(2)(p,A) 



a 
+ 



d'k9{A-f{k)) 



1 i{p — k) ■ 'y + 2m 



k^ (p — ky + 
1 i{p — k) ■ J + 2m 



J fcO = Vk2 + A2 



^0 _pO ^2 



It can be easily checked that, at variance with the expressions ( |5.12|) and (|5.14|) for the 
classical self-energy of a point charge, in the integrand of the above equation the leading 
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2 

terms of order |k| cancel each other, whereas the non- vanishing sum of the terms of order 
|k| gives rise to a contribution of order a In A to the electron self-energy. Considering 
the case < 0, we shall initially assume to choose the reference frame for which p = 0, 
= \f—p^- Taking /(k) = |k| //?, where /3 is an arbitrary parameter, the integrals can thus 
be easily evaluated. Then, by imposing the Lorentz covariance, one can generalize the result 
for arbitrary p. This provides 



/„, 2/3A 1\ 

m 3 m + ( ■ 7 + m 

\ m 2 / 




(6.5) 



where is expressible as a power series containing second and higher powers of 

(ip ■ 7 + m) and, as it can be easily checked, coincides with the finite one-loop electron 
self-energy provided by the usual renormalization theory. The form of Eq. (|6.5| ) obviously 
suggests the generalization of the notion of 7?.-limit which may be applied in Eq. ( |6.4| ): given 
a function -F(A) having a behavior for of the type 

F(A) = clnA + ci + 0(A-^) for A ^ oo, 

we shall define 

7^ lim F(A) = d. 

We obtain in this way 

yP-\p) = A(') + £(2) . ^ + ^) + (6.6) 

with 

2n \ m 2J 
5(2) = _ ln---21n — + - . 



One can see from Eq. ( |6.6| ) that A^"^"* represents a finite electromagnetic mass of the 
electron. If m is identified with the physical mass, the above result amounts to introducing 
in the hamiltonian a bare mass tuq = m— A(2)+0(a2). Besides, E^^) can be considered as the 
second-order contribution to a finite (apart from infrared divergences) fermion wave-function 
renormalization constant Z2 = 1 — -B*-^-* + O (a^). One may observe that the arbitrary finite 
parameter /3 only appears in the unobservable quantities mo and Z2: therefore it simply 
allows for a finite renormalization of the theory, once the ultraviolet divergences have been 
removed by the 7?--limit procedure. One can then conclude that with the present method 
the familiar results of perturbative QED at order a are actually recovered without the need 
of introducing divergent parameters in the hamiltonian. 

Of course, in order to really establish on a firm basis the correspondence between classical 
and quantum electrodynamics, one should be able to obtain such a result to any order in 
perturbation theory. However, an important difference between the two theories is given by 
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the fact that in the quantum case the full electron self-energy contains an infinite number 
of contributions of any order in a. When a manifestly covariant cut-off A is used (as in 
the Pauli-Villars regularization scheme) the contribution of order a" contains powers of 
ln(A/m) up to ln"(A/m). Therefore the simple logarithmic dependence which we have 
found for the lowest-order contribution in Eq. (|6.5| ) cannot be assumed to represent the 
behavior of the exact electron self-energy as a function of A. We have just seen how the 
form of this behavior (which in turn is strictly related to the large momentum behavior of the 
renormalized self-energy) affects the correct definition of the 7?.-limits that may appear in the 
equations of motion and the hamiltonian operator of QED. These considerations therefore 
seem to indicate that the problem of explicitly providing a finite and non-perturbative 
mathematical formulation of quantum electrodynamics is in any case ultimately related to 
the investigation of the short-distance behavior of the theory. 



APPENDIX A: ASYMPTOTIC EXPANSION OF THE FIELD 
NEAR A POINT CHARGE 



We want to construct a covariant asymptotic expansion for the retarded and advanced 
electromagnetic potentials near a point z{t) of the world-line of a particle with proper time 
r. Such an expansion will have the form 

<dvAet(^ + ^) ~ «dvAet)-l(^) + «dv/ret)o(^) + «dvAet)l(^) + • • • > (Al) 

where each term (^adv/ret)j('^) ^ homogeneous function of x of degree j, and contains the 
time derivatives d^z/dr^ up to the order k = j + 2. 

Clearly, the following two requirements must be satisfied: 

1. When fixing x = z + x, since the left-hand side of Eq. (ST) becomes of course indepen- 
dent of r, the series on the right-hand side must as a whole have the same property. 

2. For = XyX'^ = and x° < (resp. > 0), the well-known expression of the 
Lienard-Wiechert advanced (resp. retarded) potentials must be obtained: 



A'^.iz + x) = JLJ^ for = 0, x° <0 
A^.Jz + x) = for x^ = 0, x° > 0. 



Conversely, it is easy to see that if we manage to obtain a series satisfying the require- 
ments 1 and 2, then Eq. (|A1|) will hold for arbitrary x. Let us introduce the scalar 

y = y^x^ + {x ■ m)2 (A2) 

and the four- vector 

yi^ = x^' + xxu\^. (A3) 

We note immediately that 
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y^.W' = 0, 



(A4) 



(A5) 



If we put 



(A6) 



we get a term of degree —1 which clearly satisfies the requirement 2. Furthermore, for fixed 
X we have 



dr dr 



[x^ -2^(r)] = -M^(r), 



dy dy dx'^ ^ dy du" 1 / _ ^ ^ 

d.T r)T^ d.T r)n^ d.T ti \ ^'^ 



dr dx'^ dr du'^ dr y 
where the last equality follows from Eq, 



xxu^Xyii" 



X ■ u 



-X ■ u 



Therefore: 



u 



[x ■ u){x ■ u) — 



y 



(A7) 



In order that the property 1 be satisfied, (^adv/ret)o(3^) must give rise, when differentiated 
with respect to r, to a term of order x~^ that exactly cancels that given by Eq. ( |A7D . 
Furthermore, because of requirement 2, we must have (^adv)o (•^) ~ ^ = 0, x° < 0, 

for = 0, x° > 0. It is straightforward to verify that these conditions 



ret Jo 



and {A 
are fulfilled by 



«dvAet)o(^) = (^'^)o(x)±£«^ 



with 



u ■ X 



y 



X 



2y' 



[U ■ XjW^ 



(A8) 



(A9) 



(in formulae like ([A8| ) , the upper and lower signs of the symbols ± and =F refer to the 
advanced and retarded fields respectively). One can then iterate the procedure in order to 
obtain the higher order terms. In particular, the term of degree 1 is found to be 



(^^:dv/ret)l(a^) = T ^ {u-X-u\-x)—+iu- x)u^ + {u ■ x)u'' 



with 



(A^iix) 



3a;^ 2 f x'^ \ u ■ X 



x^ _ x^\ y.2 

y2 4^4 ) 
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From Eq. ( |A1|) one can directly obtain, using Eq. (|3J]), an analogous expansion for the 
advanced and retarded fields: 



with 



This provides 



(Fadv)-2(a;) = (Fret)-2(x) = -^u Ax= (F)_2(a;) 



(AlO) 



(Fadv)-i(a;) = (Fret)-l(x) 



2y2 



1 — — r ) (m ■ x)u A X — (u ■ x)u Ax — —u A ii 



^ (All) 



and 



with 



(A12) 



An , , 
— {F), = —{u-x) 
q 2y^ 



—u Au + {u ■ x)u A X 



X 



-^-l\[u-x) + +[u-x)[u-x) 



2y5 

1 / 

■— Ax + — 2 + — (u- x)u A u. 
2r 3y \ ?/2 ; 



u A X 



(A13) 



APPENDIX B: DIVERGENT TERMS OF THE NON-RENORMALIZED 
ELECTROMAGNETIC ENERGY-MOMENTUM 

Let us introduce, on the space at fixed t for the laboratory system, the spatial coordinates 
y„ corresponding to a reference system moving with the instantaneous velocity of the n-th 
particle: 

\ t '[\^^~ ' /-o^\ 
Yn = X - Z„ + - 1) \n- (Bl) 

V 

n 

The modulus of the vector is equal to the scalar yn given by Eq. (|2.14|) , so that we can 
write yn = \yn\- Inverting Eq. (pl|) , performing the translation x — z„ — > x and dropping 
the index n, one gets 
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from which it follows 



x-v = ^y-v, (B3) 



2 2/ \2 

X =y - (y ■ v) . 



(B4) 



Besides, from Eqs. ( [B2D and ( |A3D one obtains for x° = 

t = y' + {i-if-^v\ 



(B5) 



Defining the solid angle dVly such that y'^dQydy = d^y = 7(i^x, one obtains immediately 



y=r 



dQy 
47r 



-yiyj 



3 



Then, from Eqs. (p2|) , (|B5| ) and (pq ), it is straightforward to derive that 



dQ 



y=r 



y _ ^ij 2 

Xiyj = . 



(B6) 



(B7) 



Let us now apply the above formulae to evaluate the 1/r term of the expression ( [^.12[ ). 
From Eqs. ( |2.7|) and ( p.9|) it follows that the leading term of T'^'^ can be expressed as 



(T'"=)_4 



(F°)-2(F'^)- 



where we have used the fact that u Ax = u Ay. By using Eq. 
the corresponding leading term of P^^{t,r) is given by 



it is easy to obtain that 



in accordance with Eq. (|2.15|) . We have also 

2g2 



y'dy J dnyiT'')., 



r 2 , 



ATir 3 



2g2 



(47r)2y4 



(^°%2(n)-2 



? 2 2 



(4vr)2y 



(B8) 

(B9) 
(BIO) 



so that, using Eqs. (|27|) , (|B4D and (| 



d^x(r'^)_4^^(2/ - r) = 7"^ y'dy dnyiT'')., 



dy r dQy 
Air Jr y^ J Att 



71/^ - lyiv'yjV^ 



27. 



^ 7-^v^-- 
47rr V 3 27^ 



47rr I 3^ 67 



(Bll) 



to be compared with Eq. (|2.16|) . 
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APPENDIX C: PROOF OF THEOREM 1 



Let us consider a volume V containing a single charged particle. Omitting the particle's 
index we can write 



~dt 



y,em 



(t,r) 



d 



dt Jv 



I d^^T^^eiy -r) = f d^^doTf'^e{y-r)+ f rf^xr^%e(?/ - 
Jv Jv Jv 



- [ d^xdiT^'eiy - r) + [ d^xT^'^doOiy - r) 
Jv Jv 

+ J d^^T^-'d^eiy - r) = + J d^^T^'^d^ySiy - r) , (CI) 



where 



dV 



da niT^'{x,t) 



and we have used the equation dyT^^{x) = 0, which as a consequence of Maxwell equations 
holds for every x outside the particle's world-line. From Eq. (|2.14|) one obtains 



Oiy=Tr- = - Xi- Zi + {Xj - Zj)u-^Ui 
ax* y ^ 



(C2) 



doy 



dy_ 
dt 



dy dz^ dy du^ 
dz^ dt du^ dt 



+ _ zk)u\x, - z,)u^r^. (C3) 
oxJ y 



We can now perform the translation x — z — » x and use the same notation of App. A. Then, 
applying the definition (|A3|) with x° = 0, we can rewrite Eqs. (|C2|) and (|C3|) as 



OiV = — 

y 



(C4) 



doy = — (-y.y + Xj 



U^XhU^ 



Therefore, making in the last term of Eq. (^T|) the change of integration variables 
obtain 



— r 



(C5) 



we 



j^P'v,eJt,r) + ^^yit)=r' I y'dy J dn^T^'^d.ydiy 

= rj-^ f dfly (t^% - T^%M^7-i + T^°XjM^XftttS"0 • (C6) 

Jy=T ^ ^ 

In order to evaluate the behavior of the above expression for r 0, it is useful to expand 
the field near the particle as in Eq. ( p.29| ): 



F{x) = (F)„2(a;) + {F)^i{x) + (F)o(x) + F + 0{x) 



where (-F)-2, and (F)o are given respectively by Eqs. ( |A10| ), ( |A11|) and ( |A13|) , while 

according to Eq. ( p.30| ) F = F{0) = F{0). Using Eq. (|2.?| ) one can perform a similar 
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expansion of the integrand on the right-hand side of Eq. (|C6|). Then, since (-F)-2 and (i^)o 
are odd under space inversion, while {F)_i is even, it is easy to see that many terms of this 
expansion are spatially odd, and therefore give no contribution to the integral. Retaining 
only the terms with even parity, we obtain 



JtPv,emit, r) + ^^y{t) = {Wn-i{t)r-' + (Wnoit) + 0(r) 



with 



y=r 



(C7) 
(C8) 



-1 



y=r 



Xu 



7 VU. 



(C9) 



In Eq. ([C^), {T^'')_^ and {T^^)_^ are respectively the homogeneous terms of degree —4 
and —3 of the energy-momentum tensor, which can be simply expressed as combinations 
of (-F)-2 and (-F)-i. However, by comparing Eq. ( |C7| ) with Eqs. ( p.l5| )-( p.l6| ) we obtain 
directly 



{W 



-ir 



Airr 3 



-u 7 



Airr 
d 
It 



2 

^ • -1 

3-7 



6 dt 



(CIO) 

(Cll) 
(C12) 



Of course, Eqs. (|ClCI|) - ([CTl| ) can also be verified starting from Eq. ( |C8|) and carrying out a 
straightforward although somewhat lengthy calculation, which may only serve as a consis- 
tency check for the whole procedure. The integrations on the right-hand side of Eq. ( |C9|) 
can instead be done immediately using Eqs. ( p.9|) and (p^) . Taking into account Eq. ( P12|) 
this provides 

d 



P 

dt' 



V.era 



+ $^ = -qF^'^u.-i-y 



(C13) 



The equation of motion for the charged particle can then be obtained by postulating that, 
in accordance with Eqs. ( ^.17[ ) and ( |2.19| ), 



dt 



Using Eq. (^) this yields 

mu^" = qF^^u^ , 

which corresponds to Eq. ( p.20[ ). 

Vice versa, let us consider a volume V containing an arbitrary number of particles, all 
obeying the equation of motion ( ^.2(j| ). Then, by decomposing V as the union of a due 
number of subsets each containing a single particle, and then applying to each of them the 
equation (|C13| ), it is immediate to prove that Eq. ( |2.17|) is satisfied. 
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APPENDIX D: PROOF OF THEOREM 2 



1. The one-particle case 

Let us consider a tube T of arbitrary shape in four-dimensional space-time, containing 
the world-hne of a single particle of charge q. Let us suppose also that a given point P(t) 
of this world-line, corresponding to a value r of the proper time, is assigned the coordinates 
^/i _ ^//i _ Q inertial reference frames O and O'. We shall define tt and vr' as the 

three-dimensional spaces characterized by the equations = and = respectively. 
Clearly both spaces intersect the world-line at the point -P(f). Let us also define the vector 

where b is an arbitrary constant four-vector. From the local conservation of the electromag- 
netic energy- momentum tensor T^^", it follows immediately that 5^J^ = everywhere in 
the tube T except on the particle's world-line. For r much smaller than the diameter of the 
tube T, we have via a by-parts integration 

= j^d^xd^,J>'e{x^)e{-x'^)e{y - r) 
= - f d^^ J^e{-x'y{y -r)+ f d^^' J'^9{xy{y - r) 

JTrm JTnn' 

+ J^^ da J^'n^eixyi-x'^) - d^x^^e{xy{-x'°)S{y - r), (Dl) 



where n is the unitary four-vector normal to the element da of the tube's surface dT, y and 
yf^ are again defined according to Eqs. (|A2D -([A3|), and we have used the relation d^y = y^/y 
(note that here, at variance with the situation discussed in App. C, partial derivatives with 
respect to time are performed keeping constant the position and velocity of the particle). 
Similarly we have 



= d^x d^J^e{-xy{x'y{y - r) 
= I d^x J'^eix'yiy - r) - / rf^x' J'°0(-x°)^(?/ - r) 

JTrm JTrm' 

+ I da JX^(-x°)^(s'0) - / d^x^^e{-xy{x'y{y-r). (D2) 

JdT JT y 

By subtracting Eq. (P2|) to Eq. (|D 1| ) we obtain 

/ rf^x jOQi^y - r) - / rf^x' J"^e{y - r) 

'e^xy^-x"") - e{-xy{x'°) 

- ei-xyix"")] 6{y - r). (D3) 



da J^Un 

dT 

- 1 i^ J-^ 



y 

According to the arguments discussed in Sec. II, for r ^ the two terms on the left-hand 
side of Eq. ( P3|) can be expressed in the form 
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/ d^xj'^e{y-r) = K{P'')_^r-^ + Kn f S^T^" + 0{ 

JTnn JtHtt 



(D4) 



/ d^^' J'^e{y -r) = h'^ {P"')_^ r-^ + 6',7^ / d^^' T""" + 0{r 

JtHtt' JTnn' 

where, in accordance with Eqs. ( p.l5D and ( |2.16|) , we have defined 



(D5) 



47r V3" 6 ^ ' ' 



(D6) 
(D7) 



and similarly for the primed quantities. 

Let us now assume that the reference frame O is related to O' by a boost transformation 
bringing the particle instantaneously at rest at the proper time r. We then have v = 0, 
y = |x|, = 0, y'^ = x'^, = 7' (x° + v' ■ x), so that the last term on the right-hand side 
of Eq. (P3|) can be rewritten as 



- / d'x^ 

y 



(a;°)^(-x'°) - e{-xy{x'°)] 6{y-r) = r / dQ, 

J|x|=r 



dx^XiT^K. (D8) 



Again, it is here useful to expand the energy-momentum tensor as a series of homogeneous 
terms of increasing degree in x. By applying Eqs. ( p.7|) and with = 0, m° = 1, one 
obtains 







(47r)2|xrV2 

so that the 1/r term of the expression (03) is given by 



dVL, 



|x|=r 



dx^ Xi{T'')_J), 



dQ, 



Trr^ J\x\=r 47r 



(x-v')(x.b) 



247rr 



y-b). (D9) 



This result could also have been deduced directly from Eq. (p3D. In fact, using Eqs. (pJ4|)- 
( P^ ) and the relations b'^^u"^ = bt^u" = bg, b'^ = 7' (6° + v' ■ b), the term of order 1/r on the 
left-hand side of Eq. (|D^ ) can be written as 



1 



6 



247rr 



y-b). 



in accordance with Eq. ( p9| ). Thus Eq. ( |D9| ) is strictly related with the non-covariance of 
the divergent part (P'^)_^r~^ of the non-renormalized electromagnetic energy-momentum. 
Now, considering the zeroth order term of Eq. (p8|) , we have 



r / d^x\ 

'|x|=r J —yi-w' 



(ix°x,(r^)_3&^ = 0. 
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This result follows immediately from the antisymmetry of (T*'^)_3 under the change of vari- 
ables — s> —x'^. Therefore from Eqs. (P3|) - (P5|) we obtain 

b^n f d^^T'^" - bin f d^x'T'^''=f daJ^'n.\e{xy{-x'°)-e{-xy{x"^)].{D10) 

JTnn JTnn' J dT ^ J 

If a single charged particle is present in the whole space, one can consider the limit of 
Eq. ( |D1CI|) when the spatial dimensions of the tube T tend to infinity. Assuming as usual 
that T^^ vanishes sufficiently fast at infinity, the integral over dT on the right-hand side will 
vanish, so that Eq. (piO|) becomes 

b^p:^ - b'X'^a = f d^^T^" - bin f d^^'T""" = 0. (Dll) 

Jit Jtt' 

Since b was arbitrarily chosen, the above equation implies that P^^ transforms as a four- 
vector. One may notice that the particle's equation of motion has not been used up to now 
in this section. This means that, in the one-particle case, the Lorentz covariance of the 
electromagnetic energy-momentum four-vector defined by Eq. ( p.22|) is a consequence of the 
Maxwell equations alone, and is verified at any point of the world-line of the particle for any 
arbitrary particle trajectory. 

We would like to point out that Eq. ( piO| ) holds also in the more general case when 
in place of O one considers another arbitrary system of reference O" for which v"(r) 7^ 0. 
Assuming again that z"'^{t) = 0, and calling tt" the space with equation x"^ = 0, this means 
that we can write 



b':.n I d^x." T"^" - b'n rfVrO'^=/ daJ^n. 
Tn-TT" Jthtt' Jot 



{x"y{-x"') - e{-x"y{x'^) . 

(D12) 



To prove this, one has simply to write down the equivalent of Eq. ( PIO ) with O" in place 
of O': 

b,n I d^xT^" - b''n I rf^x" T"^"" = [ da rn. - 

JTnn JtHtt" JdT L 

Then subtracting this equation from Eq. (piO|) directly yields Eq. ( pi2|) , as it is easy to 
check. We are going to make use of Eq. (P12|) for the analysis of the many-particle case in 
the following subsection. 



2. The many-particle case 

Let us now assume that particles are present, with A^ > 1. One can prove that also 
in this case the total energy- momentum of the system, defined according to Eq. ( p. 211) , 



transforms as a four-vector when particles and fields satisfy the equations of motion (|2.20|) 



and ( |2.24| )-( p.27| ). To show this, let us divide the whole space-time into A^ tubes Ti, T2,. . .,Tj 



N 



each containing the world-line of a single particle (clearly some tubes will extend spatially to 
infinity). Given two arbitrary inertial frames O and O', the spaces vr and vr' (again identified 
by the equations x° = and = respectively) will intersect the world-line of the generic 
n-th particle at two distinct points corresponding to the values r and r' of the proper time. 
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Calling x° = x°(r'), we introduce also the space vr characterized by the equation x*^ = aP. 
Repeating the same arguments of App. C we can write 



9T„ 



(D13) 



(note that the expression in square brackets on the right-hand side is written in such a way 
that the above formula holds in both the cases of > and aP < 0). Besides, rewriting Eq. 
( P12| ) for the n-th particle at the proper time r', we have 



da J^rin 



(D14) 



Multiplying Eq. (^) by and summmg the result to Eq. (pT^) we arrive at 





-K 







JT„nTT' 



da J^n. 



Finally, let us sum over n. The integrals over the various 9T„ will cancel out (neglecting as 
usual the contributions at infinity), so that we are left with 



0, 



where the total momenta and P"^ are defined according to Eq. ( ^.21 ) in the two reference 
frames O and O' respectively. Again, due to the arbitrariness of 6, the above equation implies 
the Lorentz covariance of the total energy-momentum four-vector. 



APPENDIX E: VARIATION OF THE LAGRANGIAN 

We shall consider here a single charged particle, since the generalization of the results 
to an arbitrary number of particles is absolutely trivial. We observe preliminarily that, by 
differentiating Eq. ( p.l6| ) with respect to 5z, 6u, 6u etc., it is possible to expand the variation 
6A of the potential in a similar way: 



6A{z + x) = {6A)^2{x) + (M)_i(x) + {6A)o{x) + 6A{z) + 0{x) 



(El) 



As a function of x, {SA)j{x) is a homogeneous term of degree j which under the transfor- 
mation X — >■ —X is odd for j = —2, j = 0, and even for j = —1. Moreover SA{z) can be 
expressed as 



SAf" (z) = 7^ hm 



dQ^ 



We have in particular, using the variables defined by Eq. (|A3|) , 



(E2) 
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{5A^')-2{x) 



Airy 



(E3) 



In expressions such this, it is always assumed that Sz"^ <^ y, since in the following we shall 
first evaluate the derivatives of L(r) for fixed r, and only at a later stage we shall consider 
their behavior for r ^ in order to obtain the corresponding derivatives of the renormalized 
lagrangian L. 

Let us first consider the electromagnetic part of the lagrangian. Since the fields satisfy 
Maxwell equations, we can equivalently replace in Eq. (|3.8|) the factor ^ (|x — z| — r) with 
9{y — r), with y defined according to Eq. (|2.14| ). As it is easy to see using Eqs. ( p.8|) -( pnoD , 
Eq. (|3.1CI|) will remain valid, with only a modification in the coefficient which plays 

no role at all in our calculations. Using this new definition of Lem we have 



d'^F^''d^6AJ{y - r) - - / d'^F'^'^F.Jy S{y - r). 



(E4) 



Integrating by parts, the first term on the right-hand side can be rewritten as 

d 

+ J d^^d^F''''6A,e{y-r) 



J d'^F>^''d^6AJ{y-r) = J d'^F^^AjOiy 



+ J d^^F^''5A,d,y5{y 



Using the expansions ( p. 29 ) and (^T|) we can write for r ^ 

j d^xF°^Aj9{y -r) = ^ + nj d^xF^^Aj + 0{r 
where, according to Eqs. (p.9|), ([E3D and 



(E5) 



(E6) 



H 

r 



J d'^{F'^)_2{SA,)_,e{y-r) 



127rr 



The second term on the right-hand side of Eq. ( |E5|) vanishes because of Eq. (|3.14| ), while 
the third can be evaluated making use of Eqs. ( |C^ ) and (|C5|). After a shift x — z ^ x in 
the integration variables, followed by the transformation (|B2|), we obtain 



d^^F^'SA^d^ySiy-r) 



dQ. 



7 Jy=r 



F''^(z + x)M,(z + x)y, 



+F°\z + x)Mi(z + x) (^-VjV^ + XjV^XkU^ 



(E7) 



We can now substitute into Eq. ( [E^ ) the expansions ( p.29| ) and ([ET]), and then discard all 
those terms that have vanishing integrals due to their antisymmetry under space inversion. 
In this way we obtain 



rf^xF^^M,a^l/5(l/-r) = - + - 



dVL^ 



r 7 Jy=r 



iF"')_25A, + F'" {5A 
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where C is a term linearly dependent on 5z and 5u, which we do not need to evaluate 
explicitly. The integral on the right-hand side can be easily calculated by making use of 
Eqs. (|23), (0) and (^). This gives 



r Jy=r A-K 



- - ^ (u-'dA, + If>^'u,6z,) + 0(r). 
r 7 V 3 / 



(u'x'' - u^x') 5A^ + F^'^u^ {^Xj + Xku'^Uj) 5z^ 



ym' 



Oir) 



(E8) 



The second term on the right-hand side of Eq. ( ^4| ) can be evaluated with a very similar 
procedure, using the relation 

1 r ■ 1 ■ 1 

5y = — Xi — Zi + {xj — Zj) u^Ui bz^ H — [xj — Zj) [xi — Zi) Su^, 
which follows from Eq. (|2.13| ). We obtain 

- i I d'^F^''F,Jy5{y - r) = ^ J ^^dil^F'^'^iz + x)F^,(z + x) {y^6z^ - 
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D r 

r 27 
D q 

r 7r^ Jy=r Att 



XjU^Xidu' 



y=r 



dny{Fn^2F^^u%M + 0{r) 



da 



^u^x'F^,y^5z' + 0(r) 



- + ±F^'u^5z, + 0{r). 
r 3'y 



Therefore, combining Eqs. ([E6|) , ([E8|) and (^), we obtain 
6L,^ = 7^1im5Lem(r) = (^/ d^^F^^6Aj 



r^O ""^ ' dt 

Considering now the interaction lagrangian, we have from Eq. ( p.9| 



(E9) 



(ElO) 



(5Lint(r) = qSv' 



dQ^ 

|x|=r 4:71 



Ai (z + x,t) + g7-V 



dVLy, 

|x|=r A-K 



5A^ (z + X, t) 



dQ^ 



diA^ (z + X, t) 6z' 



'|x|=r Air 

Keeping in mind Eqs. ( p. 61) , ( p.20|) and (|E2|) we can then write 

5Lint = 7^1im5Xi„t(r) = qSv'Ai{z) + q-f~^u''5A^{z) + q-f'\''diA^{z)Sz' . 

r—*0 

We have finally 

Therefore, summing Eqs. ( ^ToD , (^) and (|ET3|) , we obtain Eq. (|02|) . 



(Ell) 



(E12) 



(E13) 
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APPENDIX F: TRANSLATIONAL INVARIANCE AND CONSERVATION LAWS 



Let us consider the Eq. (p.22|) , which has been derived making use only of the field 
equations. When also the equation of motion (|3.24|) for the particles is satisfied, the second 
expression in square brackets vanishes and we have (considering a single particle to simplify 
the notation) 



5L 



d_ 

di 



(Fl) 



Let us now examine the space translation of the whole system by an infinitesimal length 
de along the direction x'': 

dz" = 6lde, 
SA^ix) = -dkA>'{x)de. 

Substituting these expressions in Eq. ([Fl|), the manifest invariance of the lagrangian allows 
us to write 



Q_dL_d_ 
de dt 



muk 



+ qAk{z)+n J d^^F^'dkA, 



dPk 
dt 



with 



Pk = mUk + qAk{z) +11 j d^xF^WkAj . 



(F2) 



We are now going to show that the conserved quantity is equal to the fc-th component 
of the total momentum defined by Eqs. (|2.21|) and ( p.22| ). We have in fact from Eq. 

7^ J d^:>cF^^dkAj = n j d^^T"^^ + n j d^^F^^djA''. 



(F3) 



But 



J d^:siF°^djA^9{y - r) = - J d^^F^^ A^dj9{y ~ r) = - j d^^F°^ A^'^Oiy 



7 Jy=r 



y,.(F°^-)_2 [(^'')_i + A^{z)] + 0{r) 



dQy 

Jy=r An 



■47rr 



A^'iz) 



+ 0(r) 



4:7ir 



qAf'{z) + 0{r) 



so that 



7^ J d^xF^^djA^" = -gA^(z). 



(F4) 



Therefore Eq. becomes 



= mu" + n / rf^xT 
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in accordance with Eqs. ( p.21|) and ( |2.22|) . 

Let us now instead consider an infinitesimal temporal translation of the system's evolu- 
tion: 

Sz" = v'-dt, 
SAf'ix) = doA^'{x)dt. 

Since the lagrangian L contains no explicit dependence on t, we obtain from Eq. 

(mui + qAi{z)) v'-TzJ d^xF^^^doAj 



dL 
'dt 



d_ 

dt 



or 



dP^ 

IT 



with 



po = (mui + qAi{z)) -n j d^xF'^^doAj - L 



mu 



n j d^^T^^ + qA\z) +nj d^xF'^^djA^ 



(F5) 



where the last equality follows again from Eq. (0). Therefore is equal to the total energy 
that was introduced in Sec. II. In the form given by Eqs. ( [I''2| ) and ([I''5|) , the four- vector P^^ 
can also be immediately identified with the canonical energy-momentum defined in Sec. IV. 
The derivation of the conservation laws carried out in the present section clearly represents 
a generalization of Mother's theorem applied to the lagrangian defined by Eqs. (|3.2|)-(|3.5|). 



APPENDIX G: LARGE WAVE- VECTOR EXPANSION OF THE FIELD 
GENERATED BY A POINT CHARGE 



By inverting Eqs. ( ^.6|) -( pfD one obtains 

(a'^)^. (k, t)= f d^^ e-^''-^ [k' (A^)^ (x) + I (9oA'^)^._i {x)] , (Gl) 



where the particle index n has been omitted to simplify the notation. Therefore Eqs. 
( p.5| ) can be obtained in a straightforward way by substituting into Eq. ( |G1| ) the explicit 



expressions derived from Eqs. ([A6|) and (|A9|) , and then carrying out the Fourier transforma- 
tion. However, one can obtain Eqs. ( p.3D -( p75D in a shorter way by considering, for a given 
time i, the field (analogous to the "Wentzel's field" of Ref. 0) 

A'^ix) ^-J dSD{x - y)r{y)e {t - y') , (G2) 



where 

Ijt^ f+OO 

fiy) =q—S{y-z)=q dTU>^{T)5' {y - z{t)) (G3) 

7 J-oo 
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is the current density carried by the particle, while D is the usual invariant function: 



Dix) 



(2vr) 



2n 



(G4) 



It is easy to see that v4j (x) = A'^^^_{x) for [x — z (t)]^ > 0, so that we have in particular 



(x, i) = A^,, (x, i) 



(G5a) 
(G5b) 



for every x. From Eqs. (p2|)-(|G^ we get 



A^{x) = -q / d^yD{x - y) dr u'^{t)5' {y - z{t)) 



iq 



(G6) 



where r is the particle's proper time for which (r) = t. We can write for r near r 

z{t) = Z + US + ]^us^ + O(s^) 
u{t) = u + us + 0(5^), 

where on the right-hand sides z, u and u stand for z (j), u (r) and ii (r) respectively, while 
s = T — T. Substituting the two above expressions into Eq. ( |G^ ) we obtain 



A'^ix) 



iq 



(27r)3 



ds e 



-ik-us 



1 - iA; ■ M— + O(s^) 



iq 



(27r)= 



j d^k6{k^)e{k\"'<''-'^ 



+ 



„ /c ■ -it 



k ■ u {k ■ uy 



2{2tiY J k^ 



ik-{x—z) 



.{k ■ u)u^ — {k ■ u)u^ ^ ^ /^n 1-3 



k ■ u 



+ 1- 



(k-u) 



{k-uf 
Oflk 



+c>(ikr 

+ c.c. 



Keeping into account Eqs. (G5), a comparison of the above equation with Eqs. (5J.) provides 
directly Eqs. ( ^.3D -( ^^ . Of course, it is implicitly assumed that the incident field y4(^ = 
A^^ — I]n^n,ret IS Sufficiently well-behaved, so that its contribution to a'^(k, t) for large k is 



at most of order Ik 



-3 



APPENDIX H: PROOF OF THEOREM 4 



Taking into account Eqs. ( [4.7|) and (|4.16| ), it is clear that Eqs. (|5.9| ) - (|5.10|) are equivalent 
to the equation 



n j (i3xr:f^(A,7r) = -^7^ j ^W(k,t)a:(k,t) 



!(27r) 



(HI) 
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with 



1 r 



doA^doA^ + d^A^d.A^ + (A'd.A' - A^d.A") 



(H2) 



Eq. ( pi| ) can be rewritten as 



with 



r^O A— »oo 



TV 
n=l 



(H3) 



(H4) 



(H5) 



and 



1 rrf^]^ 



-k^'a''{k,t)al{k,t)e (^A - k'^ 



(H6) 



In order to evaluate the two 7^-hmits in Eq. ( P4|) , it is convenient to decompose the field 
in the following way: 



TV 



n=\ 
N 



(H7) 
(H8) 



n=l 



The two above equations provide the definition for (d^A^)^^^ and {a,'^)^es^ which are of course 
related to each other by an equation similar to ( [5.71 ): 



(duAn,^, (x) 



2(27r)3 7 P 



- jkx 



(H9) 



Let us first consider the expression ( [H6| ) . Taking into account Eq. ( [H8| ) , we have for 
A ^ cx) 



K^^^ = C^A + + 0(A 



-In 



(HIO) 



with 



2{2nfC^A =Y.J -^k^ ^ (A - k') 



and 
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n=l 



+ 



-j^k^ (a^)res (a.):es 



(Hll) 



It is easy to check that the integrals appearing in the expression of D'^ are indeed convergent. 
One has to observe in this respect that, for |k| — > oo, 



N 



(k,t) = ^«)o(k,t)e-''-- + |k 



1-3 



n=l 



and that, result of Eqs. (^-(^, 

(an,.); + K)o = 0. (H12) 

Considering now Eq. ( |H5| ) , we shall explicitly consider only the temporal component, 
since the other three can be treated in exactly the same way. We first observe that the 
three-divergence on the right-hand side of Eq. ( P^ ) gives no contribution to Il^^{t,r) for 
vanishing r. In fact, assuming as usual that the fields vanish sufficiently fast at infinity, we 
have 

N 



^ •' n=l 
n=l •' 



r). 



Each of the terms of the sum for n = 1, . . . , on the right-hand side of the above equation 
has the form 



where the index n has been dropped and a transformation of the space coordinates of the 
form (p2|) has been carried out. Using the asymptotic expansions ( |3.16| ) and (|3.19|) , with 
= qu'^/4:7iy and {d^A^)_^ = qy^u^/iny^, in the limit r — > the above integral 
becomes 



r 

'2^ 



y=r 



r 



— I dn. 

27 



qu 



{d,A^)_^-[{A^U + A^] {d,A^)_A¥ 



y=r 



+ A' 



qu-' 
Anr 



+ A^ 



Therefore, by making use of Eq. (P7|) , we obtain that the behavior of Il^^{t, r) for r — 
has the form 



r 



(H13) 
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with 



N 



r 2 ^ 



n=l 



mJ-2 



N 



n=l 



(H14) 



The fact that D'^ is a finite quantity can be as usual verified using the asymptotic expansion 
( p.l9| ) and observing that 



d'^ {dxA>^)_, e{yn - r) = 



(H15) 



because of the spatial antisymmetry of the integrand. It may be noted that Eqs. ( [H12| ) and 
( H15|) are directly related to each other. They are respectively the reasons why there is not 
a term of order Inr on the right-hand side of Eq. ( |H13| ) and a term of order In A on the 
right-hand side of Eq. ( |H10|) . 

Finally, by substituting Eqs. (13) and (^) into Eq. (|HTD , carrying out the integrations 
over the x space and comparing the result with the right-hand side of Eq. ( [Hll| ), it is 
straightforward to verify that = D'^ . Keeping in mind the definition of the 7^-limits, the 
Eq. (^^) is therefore proved. 

To complete the demonstration of the theorem, we are now going to prove Eq. ( ^.llj) . 
We first observe that, according to Eqs. ( |3.16D and ( ^.17[ ), we have 



A^' (Zr. 



lim 



A^{z^^x)-{Al)_^{x)-{Al\{x) 



so that, taking into account Eq. ( pl6|) , we get 

1 y, „ 



A^ (Zr, 



2(27r)= 



A;0 



''njR 



k,t)e 



<r^(k,t)e 



— ik-z„ 



(H16) 



with 



On the other hand, we have for A oo 
with 



eiA-k'']= E^A + F^ + 0(A-^), (H17) 
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1 



2{2Tif J fco 
1 r d^k 



(2vr) 



K)_jk,t) + «r_jk,t) 

,^j_,{k,t)e{A-k'), 



A-k' 



1 



d^k 



<)^(k,t)e^^- + K);,(k,t)e- 



(H18) 



2(27r)3y P 

One has to observe that (a(^)Q does not contribute to Eq. (pi7|) since 

«)„(k,t) + «);(k,t) = 0. 
This corresponds to the fact that (v4(:J)q is an antisymmetric function of x, and so does not 



contribute to the integral on the right-hand side of Eq. ( 3.17 ). Taking into account Eq. 
( [HlTp , the comparison of Eq. ( [H18| ) with Eq. ( |Jil6| ) provides finally Eq. ( [5.11| ). 



45 



REFERENCES 



[1] P. A. M. Dirac, Proc. R. Soc. London, Ser. A 167 (1938), 148. 

[2] F. Rohrlich, "Classical Charged Particles," Addison- Wesley, Reading, MA, 1965; F. 
Rohrlich, in "The Physicist's Conception of Nature" (J. Mehra, Ed.), Dordrecht- 
Holland, 1973. 

[3] A. O. Berut, "Electrodynamics and Classical Theory of Fields and Particles," 2nd ed., 
Dover, New York, 1980; S. Coleman, in "Electromagnetism: Paths to Research" (D. 
Teplitz, Ed.), Plenum Press, New York, 1982; P. Pearle, ibid. 

[4] A. D. Fokker, Z. Phys. 58 (1929), 386. 

[5] P. A. M. Dirac, Ann. Inst. Henry Poincare 9 (1939), 13. 

[6] F. Rohrlich, Phys. Rev. Lett. 12 (1964), 375. 

[7] C. Teitelboim, Phys. Rev. D 1 (1970), 1572; 3 (1971), 297; E. G. P. Rowe, iUd. 18 
(1978), 3639. 

[8] E. G. P. Rowe, Phys. Rev. Lett. 53 (1984), 1779. 

[9] C. Teitelboim, Phys. Rev. D 4 (1971), 345. 
[10] C. Teitelboim and C. Lopez, Lett. Nuovo Cimento 2 (1971), 225. 
[11] F. E. Low, Ann. Phys. (N.Y.) 266 (1998), 274. 



